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SUMMARY 


The K-shell photoeffect is studied in the limit of zero kinetic energy of the outgoing photo- 
electron. The relativistic (« Z)? corrections are obtained in the general case, and “exact”? nu- 
merical calculations are performed for Z=92. The total cross section is reduced by 21 % (for 
Z=92), and the angular distribution is changed somewhat, compared to the non-relativistic 
results. Considerable electron polarization effects are obtained. 


Introduction 


In the preceding paper (Nagel, 1960, hereafter denoted by I), one of the 
present authors studied the K-shell photoeffect and derived expressions for the 
matrix elements to relative order «Z in a Coulomb field, using a final state 
wave function, exact to order (x Z)*. For heavy atoms this approximation is 
poor, and one has to use the exact expression for the final state wave func- 
tion as a series of partial waves. The reaction amplitudes, in terms of which 
the photoelectric matrix elements were expressed in I, are then obtained in 
form of series expansions in radial multipole matrix elements, corresponding to 
a retardation expansion of the plane wave of the radiation field. The radial 
Coulomb matrix elements must in general be computed numerically, and as a 
large number of these matrix elements contribute, especially at higher energies, 
it is important to have suitable and rapidly convergent series expansions. This 
problem has been studied by one of us in a number of papers (Olsson, 1959; 
Olsson and Hultberg, 1959; the last paper will be denoted by II). 

At the K-shell threshold, where the photon energy is equal to the binding 
energy of the K-electron, one should expect the expansions in multipole mo- 
ments to converge rather rapidly, as the photon wavelength is then, even for 
heavy atoms, considerably larger than the radius of the K-orbit. The limiting 
form of the radial matrix elements has been derived in II in connection with 
a study of the internal conversion at low energies. In the internal conversion 
matrix elements the radiation field enters in the form of spherical Hankel func- 
tions, whereas in the photoeffect one has spherical Bessel functions. The results 
in IL can then be taken over directly, our matrix elements being the real parts 
of those in II. 
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The radial matrix elements are then expressed in terms of a confluent hyper- 
geometric function of two variables. For the computation of this function a 
rapidly convergent series in ordinary confluent hypergeometric functions can be 
used. We first derive the expressions for the reaction amplitudes to relative 
order («Z)? (the terms of order « Z vanish in this limit), and then also per- 
form a numerical calculation for the case Z=92, taking into account transi- 
tions to final states up to and including the g-states. 


1. The radial matrix elements in the zero energy limit 


The relevant radial matrix elements are defined in I (3.8). We want to ob- 
tain analytic expressions for these matrix elements for a pure Coulomb field, 
especially in the limit B=0 (B=v/c, where v is the velocity of the photo- 
electron). This is the same problem as treated in II, except that we have 
spherical Bessel functions instead of spherical Hankel functions in the integrals. 
The normalization of the radial electron wave functions is also different in I 
and II; the functions used in I (3.8) are obtained from the corresponding bound 
and continuum functions N,;% and N;Z in II (11), (2), and (24) by multiplica- 
tion with 


a 2a 
—t"' and a ee t-1, respectively. 1.1 
Vn re P y (1.1) 


Here t=2ar, a=aZ; e« is denoted by EH in II. We use in general the same 
notation as in I. 

Our matrix elements can now be expressed in ReJ, and ReJg, II (13), (14), 
and (15). We then perform the transition to the limit B=0, see II (19). The 
final result is 
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% (which is different from the x used in I) here characterizes the final state, 
%=l, or —1,—1. 

H, is a confluent hypergeometric function in two variables, and can be ex- 
pressed as a series of ,F, functions: 


aan if, (a+2m; c; y), (1.5) 


with (b)n=6(b+1) ... (6+n—1), (b)=1. 


In the expressions I (3.7) we also need the asymptotic phases 7,; for the 
Coulomb field we have (cf. Olsson, 1959, pp. 137, 138, (22) and (23)) 


p (% ea 


iaeelik (z= or —?%—1). (1.6) 


Nu=30 (l— yx) — arg D (y,+ta/p) —are tg 


The value of arctg should be chosen in the interval (0, 2). In the limit B=0 
(the arctg then goes to $2) we get, using Stirling’s formula for the I'-function, 


lim (y,.+ Mo) =3% (l—2 yx), (1.7) 
£0 


where 9) =4/B [In (a/B)-1] +42 (1.8) 


is infinite in the limit, but independent of x, so that it gives a common phase- 
factor to A, H, F, and G in I (3.7). 


2. The reaction amplitudes including order (a Z)° 


From the expressions for A, B, C, and D to first order in a, given in I 
((5.3) to (5.6)), we get, if we let 6 go to zero and define 


A’ (9) = lim i exp (i) p? A (8) (2.1) 
p->0 


and so on, gp) as in (1.8), the following limiting values: 
A’ (0)=K sin 0 exp (—ia cos 0) 
E’ (0) =F’ (0) =@' (6)=0 
with K=82V2e~ (e= basis of natural logarithms). 
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This gives the ordinary non-relativistic sin* 9 cos’ angular distribution, no 
polarizations, and a total cross section (per K-electron) 
256 2” ( h ) me 


a - = 1.679x 19-2 cm’, (2.3) 
3 et mc) hy hy 
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We will now use the expansions I (3.7) and the limiting forms (1.2) and (1.3) 
of the radial matrix elements to derive expressions for the amplitudes A’... G 
including order a®. As z is of order a, this means that we must take into ac- 
count all radial matrix elements containing j,(kr) with 1<2. The calculations 
are rather lengthy, but perfectly straightforward. For J=1 and 2 we can put. 
y,=1 and y,=|x| everywhere, both in the matrix elements and in the phase- 
factors exp (in,), but for /=0 the terms of order a? in y must be kept, and 
one then also has to expand the I-functions. We give here the resulting ex- 
pressions for the reaction amplitudes: 


A’ (0)=K [1+ (—0.1789 + 72/3) a*] sin 6 exp (—ia cos 6), 


(0) 
E’ (0) =K (0.1884 —i/12) a® sin 8, 
} (2.4) 
F’ (6)=K (0.0116 +i 2/12) a? cos 6, 
G’ (0) =0. 


K is given in (2.2). The decimal numbers in (2.4) are sums of fractional num- 
bers and terms containing Euler’s constant (from the logarithmic derivatives 
of the I’-functions). 

From (2.4) we can compute the various quantities of interest (see I (1.14) 
to (1.20)). The total cross section to order a” is 


o =0, (1—0.36 a”). (2.5) 


The relative angular distribution of photoelectrons from linearly polarized pha 
tons is 


o (6, 6) ~ sin? 6 cos* ¢ + (0.069 + 0.035 sin? 6) at. (2.6) 


(2.6) gives the correct lowest order expression for the term that is independent 
of cos’ ¢. For unpolarized photons this term should be dropped in a consistent 
approximation, and in this case we still have the non-relativistic distribution sin? @, 
The a* term gives a small emission of photoelectrons in the plane perpendicular 
to the polarization plane: 


ii 
[4m sin 9.40 / | o (0, 0) sin 646=0.14 at. (2.7) 


This term also gives: the lowest order non-zero forward and backward emission. 
intensity, as the cos’ term must always vanish for 9=0 and az (ef. I (1.14), 


(1.7), and (1.8)). We then get for the ratio of forward and backward to maxi- 
mum emission intensities 
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o (0)/o (2) =a (2)/o (bx) =0.14 a". (2.8) 
For the polarizations P,, P,, and P, (see I (1.17)-(1.20)) we obtain 


P, = —2/6a’ cot 6, 

P,=0.40 a? cos 6, 

sin? § — 0.058 
in 6 


(2.9) 
P,=0.40 a? 


3. Numerical calculation for Z = 92 


We have also made numerical calculations of the matrix elements in (1.2) 
and (1.3) for Z=92 (a=0.6715), taking into account all transitions to final 
states up to and including g7j2 and go (x=4 and —5). The H, functions were 
calculated from (1.5), where the ,F, functions are obtained from their power 
series expansions. 

The resulting radial multipole moments 


a, = lim p' a,, (3.1) 
p—0 


and similar for bz, are given in Table 1. From this table we see that the 
magnitude of the multipole moments decreases rapidly with increasing multipole 
order; the #1 multipoles are dominant, as could be expected, and next in im- 
portance comes #2, and then M2 and M1. 


Table 1. Table of radial matrix elements a, and b, (3.1), corresponding to ab- 
sorption of electric and magnetic multipoles under transition from the initial 
1s, state to final state, specified by x. Z=92, B=0. 


Transition to final Electric , Magnetic , 
state with x multipole 4g multipole = 
$1, =i) “=: = M1 — 0.0106 
Pr), 1 Fil 0.0963 —s — 
Ps), —2 El — 0.2283 M2 0.0136 
ds), 2 E2 0.0286 M1 ~0.0017 
ds), -3 B2 — 0.0390 M3 0.0034 
fs), 3 E3 0.0038 M2 -0,0002 
fry, —4 E3 — 0.0046 M4 0.0005 
7), 4 B4 0.0004 M3 0.0000 

; —5 BA — 0.0005 M5 0.0000 


Using the results in Table 1 and the phases obtained from (1.7) we have 
computed the functions A’ (6) ... @’ (8), and from this tabulation the cross see- 
tions and polarizations can be calculated. 

The total cross section is obtained directly from Table 1 and I (3.9), and is 


= 0789.65: (3.2) 
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Fig. 1. Angular distribution of zero-energy K-shell photoelectrons from Z=92. The curve 
drawn with thin line is the non-relativistic sin? 9 distribution. The curves are normalized to one 
at 0=90°. 


Fig. 2. Polarizations of zero-energy K-shell photoelectrons, Z=92. P, is a transverse polariza- 

tion perpendicular to the scattering plane, and is obtained also for unpolarized photons, P, and 

P, are longitudinal and transverse polarizations in the scattering plane, corresponding to circu- 
larly polarized photons. 


6, is the non-relativistic value given in (2.3). We thus get a reduction of the 
cross section with 21% (the a* term in (2.5) gives 16%), provided we use in 
both cases the same value for the photon energy hv. The factor mc*/hy comes 
from the normalization of the incident photon flux, so it is not entirely in- 
consistent to use both relativistically and non-relativistically the same (experi- 
mental) value, although the theoretical binding energies in a pure Coulomb field 
are rather different in the Schrodinger and Dirac cases (115 and 132 keV, re- 
spectively; the experimental value is 116). 

The angular distribution of photoelectrons from unpolarized photons is shown 
in Fig. 1. The distribution is pushed somewhat towards larger angles (cf. dis- 
cussion in I, Sect. 6a), the maximum emission occurring for 6=97°. About 
45 % of the photoelectrons are emitted into the forward hemisphere (@ < 90°). 
We also get finite forward and backward intensities, amounting to 1 and 7 %, 
respectively, of the maximum intensity. The asymmetry around 90° is mainly 
due to interference between the outgoing p- and d-waves. 

If the photons are linearly polarized, we get the value 0.05 for the ratio 
defined in (2.7). 

The electron polarizations are shown in Fig. 2. The degrees of polarization 
can be rather large, especially for angles near 9=0 and 180°. For large back- 
ward angles (around 160°) one obtains, even for unpolarized photons, a trans-— 
verse polarization of more than 50%. The lowest order expressions (2.9) give 


the general trend of the §-dependence, but the quantitative agreement is in 
general not very good. 
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The relative inaccuracy, due mainly to neglect of higher order multipoles, in 
the numerical results underlying Figs. 1 and 2 is estimated not to exceed a 
few percent for the polarizations. The angular distribution is probably more 
accurate, and for the total cross section the contribution from multipoles with 
[I> 2 is less than 0.3 %. 


The results of this paper should be applied with some care to the actual photo- 
effect at the K-shell threshold. When the outgoing photoelectron has zero energy 
one could expect the influence of screening, particularly on the polarizations, 
to be considerable. The simple screening correction to the total cross section 
(see e.g. Bethe and Salpeter, 1957, p. 391; this correction assumes that one has 
an analytical expression for the total cross section near the threshold), which 
consists in replacing the actual potential in the K-shell region by a displaced 
Coulomb potential with an effective charge, is hardly applicable to angular dis- 
tributions and polarizations, as these effects are sensitive to the relative phases 
of the different partial waves. 
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